cm UK Patent Application „„GB ,,,,2 309 562 „ 3 ,A 



(43) Date of A Publication 30.07.1997 



<21) Application No 9701296.7 
(22) Date of Filing 22.01.1997 



(30) Priority Data 
<31> 9601461 



(32) 26.01.1996 (33) FR 



(71) Applicant(s) 

Instrtut Francais du Petrole 

(Incorporated in France) 

4 Avenue De Bois Preau, 92852 Rueil Matmaison, 
Cedex, France 

(72) Inventor(s) 

Philippe Joseph 
Didier Granjeon 

(74) Agent and/or Address for Service 
Frtzpatrkks 

Cardinal Court 23 Thomas More Street LONDON, 
El 9YY, United Kingdom 



(51) INT CL 6 

G06F 17/00 , G01V 1/28 , G06F 17/13 

(52) UK CL (Edition O) 

G4A AUA 

(56) Documents Cited 

EP 0297737 A2 US 4340934 A 

Inspec Abstract Accession No 02723383 & Journ. 

Geodynamics, Vol 3, No 1-2, July 1985, pp 155-168 

Inspec Abstract Accession No 03279775 & Earth 

Surface Processes & Landforms. v13,n6,Sep 88, 

P487-500 

(58) Field of Search 

UK CL (Edition O ) G4A AUA AUXX 

INT CL 6 G01V 1/28 1/30 1/40 3/38 9/00 11/00 . G06l= 

17/0017/13 17/50 

Online: COMPUTER, INSPEC, Dialog/GEOPMYS, WPI 



(54) Modelling the filling of sedimentary basins 



(57) The method relates to the construction of a 
numerical deterministic model of the diffusive type 
allowing the filling of sedimentary basins by siliciclastic 
and carbonate sediments to be simulated in 2-D or 3-D. On 
the basis of known data pertaining to the architecture of a 
basin and measurement data, well loggings, seismic 
surveys, etc., a set of input data is produced relating to the 
accommodation structure made available by subsidence 
and eustasy, the inflow and production of fluviatile or 
marine sediments and physical transportation parameters 
such as the diffusion coefficients of the different 
lithologies. This set of data is applied to a numerical 
model. The results that can be derived therefrom as 
regards the geometry and lithologies of the sedimentary 
units are compared with measurement data and the input 
data are iteratively fined down, step by step, by a process 
of inversion. The model constructed in this manner allows 
the volumetric partitioning, distortion of the genetic units 
or sedimentary geometries observed, for example, to be 
reproduced. 
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1 

METHOD FOR SIMULATING THE FILLING 

OF A SEDIMENTARY BASIN 

The objective of the present invention is a method 
for building a model to simulate the sedimentary filling 
5 of basins on large scales in terms of time and space. 

The method of the invention relates more specifically 
to the creation of a numerical strat igraphic model that 
will allow mult i - lit ho logical filling of a basin to 
simulated in two or three dimensions (2-D or 3-D), the 
10 purpose being to simulate the st ratigraphic response of 
the sedimentary system to variations over time in the 
eustasy, subsidence, sedimentary deposits as well as the 
physical parameters governing the inflow of sediments to 
the basin. 

15 BACKGROUND OF THE INVENTION 

Recent progress in the field of geology, from which 
seismic stratigraphy followed by genetic stratigraphy has 
grown, has profoundly changed our understanding of how 
basins have filled up with sediment on large scales of 
20 time and space, by demonstrating the primordial effects of 
three main parameters: eustasy, tectonics and sedimentary 
flux. 

Numerous models, in particular numerical 

deterministic models, have been created as an aid to 

25 understanding the geometric and lithological implications 
of these new approaches. 

These numerical models simulate the transportation 
and depositing (or erosion) of sediments in the basin by 
applying a more or less complex description of the nature 

30 of subsidence and the inflow of sediment based on an 
estimated eustasy. It should be recalled that the term 
eustasy is used to refer to variations recorded in the 
surface of the oceans simultaneously across the whole of* 
the earth's surface whilst the term subsidence is used to 

35 refer to the absolute displacement of the bottom of a 
sedimentary basin relative to a fixed marker level over 
time. These models therefore enable the effects of 
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different theories (the relationship between climate and 
eustasy, etc.) on the layout of sedimentary units to be 
tested. In addition, when applied to real situations, they 
enable the consistency of parameters incorporated in the 
model, such as eustatic and tectonic variations, to be 
tested and the geological interpretation of the basin 
under study to be confirmed. 

A model is expected to be able to define the major 
trends in variations of the facies (variation in the 
sand/clay ratio, carbonate content, etc J within genetic 
units. A mult i - lithological model -building system is 
needed for this purpose. The model must be capable of 
simulating the transportation and sedimentation of 
different siliciclast ic (sand, clay, etc.) and carbonate 
lithologies (reef, deep-sea carbonate mud, bioclasts, 
etc.) so that these trends in facies variations are 
generated by the model, irrespective of the given 
geological conditions, such as the sandy alluvial plain 
and clayey marine area. In this ideal model, the 
siliciclast ic sediments are brought into the basin by 
rivers and sea currents. The carbonates are produced in 
the marine area in accordance with a law based on the 
bathymetry and turbidity of the water. 
P = P Q A turb e" k2 

P productivity of carbonates 
P 0 optimum productivity 

A turb attenuation factor depending on the turbidity 

(turb.) of the water 
e" Kz attenuation factor depending on the bathymetry 

(z) . 

Once they have been brought to the level of the basin 
boundaries or produced in the basin, the sediments are 
transported and settle. 

Three major families of deterministic models govern 
the transportation of sediment: 

particulate models based on solving the motion of 

particles (calculation of the water flow followed by 
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the relationship between the water flux and the flux 
of sediments) ; 

diffusive models based on an equation of diffusion in 
which the diffusion coefficient can be more or less 
precisely defined (taking account of several 
lithologies, the water flux, the environment in which 
deposits are being made, etc.); and 

geometric models based on a definition of the 
geometric profile of the environment in which the 
deposits are occurring (length and incline of the 
alluvial plain, etc.) or on a geometric definition of 
the rate of sedimentation (exponential decrease in 
the marine area, etc.) . 

Particulate models use a careful definition of the 
sedimentary process and are therefore as chaotic as 
nature. Basically, they allow simulations to be produced 
on the scale of a reservoir (length approximately 0.5 to 
50 km and duration of 5 to 500 ka) . Diffusive and 
geometric models both provide more or less rough 
approximations of the natural processes. They are much 
more stable than the real situation they are intended to 
portray and provide only a toned-down approximation of 
nature. These models are applied by preference on the 
scale of the basin (length approximately 10 to 1,000 km 
and duration of 0.1 to 10 Ma) . 

Geometric models, based on an approximation of fhe 
morphology of basins, are mainly applied to straight- 
forward 2-D cases, in which subsidence, the nature of the 
sediments and climate do not affect the definition of the 
equilibrium profile too extensively. 

More general in nature, diffusive models, based on an 
approximation of the physics of sediments, are applied in 
3-D and are capable of processing transportation and 
sedimentation of multiple lithologies. 

Various known diffusive models are described, for 
example, by: 

- Kenyon and Turcotte, 1985, Morphology of a delta 
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prograding by bulk sediment transport, in Geol . Soc . 
Amer. Bull., 96, 1457-1465. 

- Begin, Z.B., 1988, Application of a diffusion - erosion 
model to alluvial channels which degrade due to base- 
level lowering. Earth Surface Processes and Landforms, 
13, 487-500. 

- Rivenaes, J.C., 1988, Application of a dual-lithology, 
depth dependent diffusion equation in strat igraphic 
simulation, Basin Research, 4, 133-146. 

One advantage of diffusive models is that they allow 
a return to geological concepts because they quantify 
relations such as the variable duration of progradation 
and retrogradation phases of genetic units or changes in 
the sandiness either on the basis of the bathymetry or the 
prograding or retrograding tendency of the littoral. 

The purpose of the known models is essentially that 
of numerical experimentation with theoretical concepts 
along 2--D profiles. 

The modelling method of the invention makes it 
possible to build a deterministic model of the diffusive 
type in 2-D or 3-D, which simulates the multi -lithological 
filling of a sedimentary basin. It uses known field data 
pertaining to the basin architecture and measurement data 
such as well- logging or seismic data to make up a set of 
input data relating to a space made available by 
subsidence and eustasy, the depositing of fluviatile or 
marine sediments and the transportation thereof, and 
physical parameters such as the diffusion coefficients of 
the various lithologies. The method is characterised in 
that it comprises, by a process of iteration: 

- dividing the basin into grids of regular dimensions, 

- solving large-scale diffusive transport equations in 
accordance with an explicit pattern of finite volumes 
and constant time intervals, so as to simulate the flux 
of each lithology deposited on each grid; and 

- comparing the simulation results with the field data and 
modifying the input data step by step by a process of 
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inversion . 

The inversion step may be performed by a process of 
trial and error. 

The diffusive model produced provides an 
5 understanding in 3-D as to how a basin is subjected to 
multi-lithological filling and therefore tests the 
consistency of a data base or even helps to predict the 
position and the geometry of sedimentary bodies. It gives 
a realistic prediction of the thicknesses of sedimentary 
10 bodies and facies. These parameters, adjusted to the well 
data, therefore make it possible to refine what is known 
about the basin overall by pinpointing the location and 
quality of bodies that are likely to be hydrocarbon 
reservoirs . 

15 Since, in order to re-construct the sedimentary 

geometries, the model of the invention requires that the 
physical parameters be accurately re-constructed on a 
large scale (flux quantification, diffusion coefficients, 
etc.), it is necessary to identify which links will have 

20 to be defined between accommodation structure, sedimentary 
deposit and diffusion coefficient in order to re-construct 
the geometries observed and hence the proposed pattern of 
correlation of the well data. The geological consistency 
or non- feasibility of these various links will therefore 

25 make it possible to confirm or invalidate the selection of 
one pattern of correlation as opposed to others. 

Other features and advantages of the method of the 
invention will become clear from the following description 
and with reference to the appended drawings, in which: 

3 0 - Fig. 1 illustrates the modelling process; and 

- Fig. 2 illustrates a model optimisation pattern by a 
process of trial and error. 

The deterministic model is constructed in a series of 
successive stages. In a first stage, the accommodation 

35 structure for the sediment is defined. The sediments are 
then introduced and sedimentary production in the basin is 
simulated. These sediments are then distributed in 
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accordance with the laws of macro- transport , which will be 
defined in the description below, limiting erosion to the 
weathering layer. 

1) Creating the accommodation structure 

The space available for the sedimentary filling of 
the basin is the sum of the eustasy and subsidence. It was 
decided that this would be defined point by point, either 
on the basis of eustatic curves and subsidence maps or 
directly from maps of accommodation structures, without 
therefore using any physical models linking the eustasy 
with climate or subsidence with tectonic, isostatic or 
thermal processes . 

Although it does not modify the accommodation 
structure, compaction affects sedimentary filling by 
altering the thickness of sedimentary layers. In order to 
take account of mechanical compaction, it was decided that 
the porosity of the sediments would be directly linked to 
the maximum burial reached thereby, using a relationship 
of exponential form which would provide a good 
approximation of the compaction, as defined by 
- Beaudoin et al , Mesure directe de la compaction dans les 

sediments, In: Berthon, J.L., Burollet P.F., and 

Legrand, Ph. (eds) . Genese et evolution des bassins 

sedimentaires . Notes et Memoires No. 21, Total 

Compagnie Frangaise des Petroles. 

In order to take account of the case of multi- 
lithology sediments, it was assumed that each lithology 
had a porosity independent of the others, which amounts to 
likening the sedimentary layers comprising a mixture of 
several lithologies such as sand and clay to a 
superposition of multiple sub- layers made up of pure 
lithologies. The porosity linked to each lithology is 
therefore dealt with individually as a function of the 
maximum burial reached by the sedimentary layer. 

for each lithology i =p <t> x = d> r , + ^o.i " ^r.J e* z/z . 
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where 



porosity of the lithology {in %) 



residual porosity of the lithology i (in %> 



deposit porosity of the lithology i (in %) 



I z 



maximum burial reached by sedimentary layer studied (in m) 



reference burial of the lithology (in m) . 
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In particular, this definition of porosity linked to 
the lithologies allows the individual transport of each 
lithology to be simulated and its effects on the porosity 
of the sedimentary layers to be deduced therefrom. 

II) Introducing and producing sediments 

Having created the accommodation structure in this 
way, the second stage of building the required model is to 
introduce sediments into the basin or produce them within 
the marine area. 

1) Introducing sediments into the area studied 
In geological terms, the flux of sediments at the 
boundaries of the area under study represent the deposits 
of sediment responsible for filling the basin. They can be 
physically perceived in terms of sediment flux boundary 
conditions . 

These boundary conditions are fixed either by 
imposing the exact value of the flux in one sector of the 
boundary or by imposing continuous change on the sediment 
flux. The first of these instances represents a deposit 
zone imposed by conditions outside the basin, such as the 
outlet of a river draining a river basin external to the 
basin. The second of these cases represents a free zone 
along which the sediment flow is governed by the physical 
parameters within the basin, such as the characteristics 
of the waves. When simulating a basin, it is possible to 
combine these two types of boundary conditions, making a 
distinction between a continental zone, for example, where 
the flux is imposed by external fluvial deposits and a 
marine zone in which the flux is defined by the laws of 
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internal transport . 

b) Sedimentary production inside the area studied 
The sediments may also be produced inside the basin 
and more particularly so in the case of carbonate 
5 sediments. For the purposes of this modelling phase of the 
invention, it was decided that this process of production 
would be simulated with the aid of an empirical formula 
linking the production rate at each point of the basin to 
the bathymetry and the sediment flux as defined, for 
10 example, by 

- Lawrence et al . , 1990, St ratigraphic simulation of 
sedimentary basins: concepts and calibration in Amer. 
Assoc, Petrol, Geol . Bull., 74, 3, 2 73-295. 
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In order to take account of the erosion of reefs when 
these are emerged, the process of mechanical weathering of 
these reefs was simulated by assuming that any reef 
sediment located in the weathering layer in the 
continental area would be transformed into a " bioclast ic " 
sediment. The reef and bioclastic lithologies are 
transported into the basin in the same way as the 
siliciclastic sediments . 



r 
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III) Sediment transport 

Having defined the accommodation structure and 
quantity of sediment that must fill the basin, the model 
distributes the sediments across the entire basin in 
5 accordance with the laws governing large-scale transport, 
established by building on known transport laws, allowing 
an estimate to be produced for the flux of each lithology 
flowing into the basin at all points. 



, oh 

10 I diffusion equation Q = — K ~~~ 

i 

i — -it 

L conservation of mass _ — — ^ 



15 where f Q sediment flux (in mVyear) 

i K diffusion coefficient (in mVyear) 

I h altitude of the point studied (in m) . 



The rate of sedimentation or erosion is defined by 

20 linking these laws to the principles of conservation of 
mass. However, this sedimentary transport is limited to 
the weathering layer and the erosion velocity at a given 
point may not exceed the local weathering velocity. 
Creating a large-scale transport law in 3-D 

25 The law of small-scale continental transport given 

above represents flows of water and sediments on time- 
scales ranging from a few hours to a few months and space 
scales in the order of several square kilometres. In order 
to be able to apply these laws to large scales of time and 

30 space, it was necessary to convert them into large-scale 
transport laws. 

To move from the small-scale continental area in two 
dimensions to the basin as a whole on a large scale and in 
three dimensions, it is necessary to: 

35 - take account of the diversity of flows in the 



continental area ranging from the minimum low-water 
patterns where the water is directed into channels 
through to the flood patterns where the water spreads 
across all the flood plains; 
5 - in the marine area, extrapolate a transport law 
validated for the continental area and define how the 
action of the waves will be simulated; 
- describe the differential but simultaneous transport of 
several lithologies; and 
10 - take account of the limited availability of sediment. 
Expressing the law of fluvial transport in 3-D 
In order to express this law of transport in three 
dimensions and take account of the various water flow 
patterns, from flooding down to low-water levels, it was 
15 assumed that the flux of sediment at any point was a 
continent proportional to an intrinsic diffusion 
coefficient, the flux of water and the incline of the 
basin . 

The water flux, constant during the study of small - 
20 scale transportation of sediments in two dimensions, now 
incorporates the variability over time of the different 
water flow patterns. Three basic flow models have been 
defined . 

a) The straight flow corresponds to the situation 
25 where the water arriving at a given point of the basin is 

entirely re-distributed in the direction of the largest 
incline and where the time interval separating two 
simulation instants is short relative to the period of 
over which changes occur in the course of the water. This 

30 period is equivalent to the meandering period if the basin 
is studied using extremely fine grids of a size of about 
one hundred metres and is equivalent to the avulsion 
period if the grids are about one kilometre. A water flow 
of this type gives rise to a flow of sediment in a single 

35 direction at each instant. 

b) What is referred to as "wild" flow corresponds 
more specifically to the situation of large-scale 
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transportation of sediment, where the simulation time 
period used is about a millennium. It is therefore longer 
than the period during which changes occur in the course 
of the water. For this type of flow, it can therefore no 
5 longer be assumed that an individual stream has kept to a 
single direction over the period and it becomes necessary 
to take account of the different courses that the water 
might have followed. In addition, unlike the laws 
governing small-scale transport, which focus on a single 

10 flow pattern, such as the flood pattern or low-water mark 
pattern, the large-scale transport law must take account 
of each of these regimes in order to reproduce accurately 
both the flood process, where clays are deposited on the 
flood plains, and the low-water level process, where 

15 fluviatile incision occurs and the sediments are rapidly 
transported towards the delta area. 

In the case of a wild flow of this type, it is 
assumed that the water has followed all those directions 
in space having a positive incline. However, the choice of 

20 these directions is made at each instant by taking the 
path of greatest incline. This being the case, it must be 
assumed that the probability of choosing one direction 
rather than another is equal to the ratio of the incline 
of this direction to the sum of all the inclines. The 

2 5 greater the incline of one direction relative to the other 

directions, the greater chance there is that the water 
will flow in this direction. This wild flow is the average 
of all the straight flows occurring over the course of the 
simulation period. 
30 c) Uniform flow corresponds more specifically to the 

case where the reasoning outlined above is generalised by 
applying a mega-time scale of one hundred thousand years. 
In this case, it can be assumed that the relief has 
reached a state of equilibrium. The flow of water 

3 5 therefore becomes uniform and there are no areas of more 

prevalent flow such as an incised valley or a belt of 
meandering. Water flow of this type gives rise to a 
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uniform diffusion of sediment. 

These three flow types, all combinations of which are 
possible, therefore represent three ways of looking at the 
drainage system of a basin over different lengths of time. 
These flows correspond to the changeover to the average of 
the straight flows, giving rise to wild flow, and the wild 
flow giving rise to uniform flow. 

In order to take account of the fact that within the 
wild and uniform flow modes the flux of water does not 
take a single direction but an infinity of probable 
directions, it was decided to replace the term water flux, 
which implies direction, with the term water efficiency. 
This water efficiency was then standardised by dividing 
the water flux by a reference flux. 

Extrapolation to the marine area 

A law of transport applicable to the marine area was 
also chosen. To this end, the various types of sediment 
transport were classified as a function of the origin of 
their motive energy, making a distinction between 
"gravity 11 flows drawing their kinetic energy from a 
transformation of the potential energy of the water flow 
and flows induced by the movement of waves due to the 
action of winds. 

It was verified that gravity flows of sediment could 
be simulated by a diffusion equation of the same type as 
that defined by Kenyon and Turcotte, 1986, mentioned 
above . 

Similarly, three types of flow were considered, the 
first being straight flow where the water is channelled at 
each instant, wild flow where the average water flow is 
considered by taking account of the distribution of water 
in the calculation and uniform flow which gives rise to a 
uniform distribution of sediments leading to regression of 
the littoral in the sector located near to the mouth of 
the delta. 

The choice of flow type depends on the value of the 
interval separating two simulation instants and the 



fluvial hydro -dynamism of the basin. If the time interval 
is short or the hydro -dynamism high, it is necessary to 
simulate a straight flow. A wild, even uniform flow, on 
the other hand will enable the envelope of sedimentary 
5 bodies to be represented. 

In order to construct the model of the invention, 
empirical extrapolation of the diffusion equation to the 
basin as a whole is assumed to be valid for "gravity" 
flows . 

10 Wave -induced sediment transportation 

As with fluviatile flows, it is possible to combine 
the laws governing fluid mechanics in order to obtain a 
small time-scale law (a few days to a few months) that 
will describe wave-induced sediment transport. On the 

15 other hand, it is not possible to extrapolate such a law 
on large time scales. It is extremely difficult to re- 
construct the exact morphology of the oceans at each stage 
of sedimentary filling and thus define a complete 
geological model. For this reason, empirical laws 

20 describing changes in the littoral on a large time scale 
were used, linking the sediment flow induced by waves to 
local parameters such as the incline of the basin and the 
local characteristics of the waves. 

Wave- induced marine transportation orthogonal to the 

25 littoral 

In the model, sediment transportation to the 
shoreline is simulated by a diffusion equation and is 
incorporated in the fluvial transport system described 
above with the aid of a negative diffusion coefficient 

30 representing the braking effect of the waves on fluvial 
transportation. Because of the problems of physical and 
numerical stability of the sedimentary system, it was 
assumed that the waves slow down the movement of sediment 
towards the open sea but are unable to generate a residual 

35 flux of sediment towards the coast on the millennium 
scale, which means that this braking effect can never 
reverse the direction in which the sediment is 
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transported. In order to account for the reduced activity 
of the waves in the direction of the open sea, it was 
assumed that the coefficient decreases exponentially with 
the bathymetry. 
5 Wave -induced longshore drift 

An empirical relationship linking the sediment flux 
flowing at a given point to the total energy of the wave 
motion and to an attenuation factor linked to the 
bathymetry was used as a means of quantifying longshore 
0 drift so as to use only local parameters. 

Simultaneous transport of several litholoaies 
The large-scale laws of transport mentioned above 
were established for a soil comprising a single lithology. 
In order to take account of a situation where multiple 
5 lithologies are transported simultaneously, the following 
relationships were used, incorporating relative diffusion 
coefficients : 



0 



5 



for each lithology i 




diffusion => 








flux of the lithology 


I (in mVs) 


1 lv s 


lithology I content (dimensionless) 


where < K, 


diffusion coefficient 


relative to the lithology I (in m*/s) 




water efficiency (dimensionless) 


th 


altitude of the ground surface (in m) . 




and = K Muv j 





These large scale transport equations are linked to 
the theory of conservation of mass, applied to each 
lithology so as to define at each point of the basin the 
rate of erosion or sedimentation and the content in each 
of the lithologies. 

Sedimen t availability and the weathering layer 
In order to take account of the uneven diffusion 
capability of the sediments of each lithology, it was 
assumed that, as already described by Rivenaes, J.C., 



r 
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1988m, op. ext., there is a layer of constant thickness 
(decimetric, for example) in which sediment transport 
takes place. The sediment flux is then proportional to the 
lithological contents of the weathering layer. 
5 In order to construct the model of the invention, it 

was decided that the weathering velocity, or the speed at 
which this superficial weathering layer is created at each 
point of the basin, would be estimated by means of an 
empirical relation . 
10 In view of the above, it was decided that for the 

purposes of building the model of the invention, the 
following two empirical large-scale transport laws would 
be used, the first relating to diffusive transport and the 
second to shore-line transport: 

15 



20 



25 



for each lithoiogy i 



littoral drift 



Q= Kr^Cji^sinacosar 



lb/bp if b<bp 

where f <b) * le " ,t k P' if b>bp 

|Q S flux of the lithoiogy i induced by littoral drift (in nr/s) 

lie , efficiency of the transport of the lithoiogy i by the waves 



density of the wave energy 



* rate of propagation of the wave energy 

lf(b) flux distribution function as a function of the bathymetry 
la angle of approach of the waves 



3 0 Numerical solution 

In order to construct the model of the invention, the 
numerical solution to transport equations based on a 
discrete treatment of the basin under study in space and a 
discrete treatment of the formation in time has to be 
3 5 found. 

The basin studied is divided into square grids of a 
constant size, whilst the filling of each grid is 
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simulated on the basis of a succession of calculating 
times, separated by a constant time interval. The width of 
the grids is from 1 to 10 kilometres, depending on the 
size of the basin being simulated, whilst the time 
interval is in the order of 50,000 years. 

The transport equations are then solved using an 
explicit numerical pattern in which the physical 
parameters such as flux or altitude at the time (n+1) are 
expressed as a function of parameters measured at the time 
(n) , the values of which are therefore known. 

In the case of a single lithology in two dimensions, 
this solution is expressed as follows: 

(1) G>=-A— => in each grid i \J = - K . ; 



* cq nr-tir QT"-sC 



(2) ™ =— _ =^ in each grid i f 

a ck at dx 

-e. h. =h JC, — 5T"J 

r / (,,) 

\ fit altitude of the grid i and at the time n {in m) 

\ K^ t diffusion coefficient at the level of the grid i {in m : /s) 

1 \y flux of sediments between the grids i and 1 + 1 (in m*/s) 

I dt time interval (in s) 

U arid width {in m) 



However, this rather simple form of expression 
becomes more complex as we move to a situation of multiple 
lithologies in 3-D and once the restriction on the 
transportation of sediment to the weathered layer is taken 
into account. In such a case, it can be said that: 



r 
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~ = - y t max , i.e. 

tir u =tir-v,.d< 

Explicit solution by finite volumes is the fastest 
5 calculating method. Furthermore, it provides very accurate 
results as long as it remains within the stability range 
of the algorithm. which requires very short internal 
calculating time intervals in the order of a century to a 
millennium . 

10 By applying the model to the input data set up from 

the field data, it is possible to simulate the transport 
of sediment across the entire basin being studied. 

The validity of the simulation is then tested by 
matching the results provided by the model with the data 
15 compiled in the field, mainly the thicknesses of the 
sediments and the facies observed on well logs. 

If there is a discrepancy, the set of input 
parameters for the model is sought by a process of 
inversion - accommodation structure and quantity of 
20 sediment transported - so that the difference between the 
results produced by these parameters and the constraints 
imposed are minimal, as illustrated below, 
direct model => S=M(p) 
difference function => E(p) = lls-dl 

2 5 ( p whole of the input parameters of the model 

I M whole of the equations governing the model 

{ S whole of the results at the model output 

! C whole of the geological constraints 

I E difference function between the results and the constraints 

3 0 1 llxll measurement of the different between results and constraints 

inversion =s> find p such that 

E=E(p) is the absolute minimum of the function E 
Vp, E(p)>E 

Within the context of the present method, the purpose 
35 of the inversion process is to quantify in particular the 
values of the accommodation structure, the sedimentary 
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inflow and the diffusion coefficients in order to produce 
a simulation in which the geometries of the sedimentary 
bodies as well as the measured thicknesses and facies 
directly below the calibration wells are as close . as 
5 possible to the geological constraints. 

A method of the trial and error type is used. A set 
of initial parameters is defined. The solution to the 
model linked to these parameters is calculated and the 
value of the parameters modified as a function of the 
0 discrepancy between the solution and the geological 
constraints. This process is continued until the 
discrepancy becomes sufficiently low. 

Once the optimum concordance has been reached, the 
simulation will provide quantitative data pertaining to 
5 the geometry of the basin and the lithologies of the 
sedimentary units. It will also provide a means of 
checking that the correlation pattern of the wells is 
consistent . 

Geological validation of the model obtained; 

The validity of the model produced in this manner was 
tested on simple two- and three-dimensional theoretical 
examples. It was assumed that the subsidence velocity, the 
rate of sedimentary inflow and the diffusive action of the 
sediments were constant. By defining the eustasy in 
accordance with one or two orders of cyclicity, these 
simple but realistic examples were used to demonstrate 
that the diffusive model is capable of representing the 
implications of theories of genetic stratigraphy, such as 
the volumetric partitioning of the sediments within a 
genetic unit and the distortion of the genetic units. The 
geological consistency of the results of these simulations 
is such that it provides empirical validation for the use 
of the multiple- lithology diffusion equation on large 
time and space scales. 

The model was also tested on real situations. The 
size of the sedimentary objects simulated was typically on 
the scale of a reservoir - Mesa Verde Group of the San 
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Juan basin (Colorado, USA) - for example, and on the scale 
of a basin - the lower Cretaceous of the Paris Basin - and 
the inversion method was used to define the parameters of 
the model (accommodation structure, sedimentary inflows 
5 and diffusion coefficients) based on geological 
constraints (thicknesses and facies of the sedimentary 
layers) by means of the trial and error method. 

This inversion, coupled with the non-dimensioning 
theory of the diffusion equation produces simulations from 
10 which the results are adjusted to the geological 
constraints. It was verified that the diffusive model used 
does produce very accurate simulations, with mean 
differences for the cumulative thicknesses of the 
sediments of some 5 metres to 25 kilometres away from the 
15 constraint wells (for a formation of an average thickness 
of 100 metres) and a difference of less than 10 kilometres 
(the size of the calculating grids was 10 kilometres) 
between the position of simulated and observed shores. 

As a result of this geological and faciological re- 
20 construction and the fact of having access to the physics 
of sedimentary processes on a large scale, the diffusive 
model of the invention allows a return to the geological 
data bases in that it fines down the values of the 
bathymetry and sediment flux, for example, by re-producing 
25 changes in the accommodation structure over time and by 
confirming or invalidating the selection of a pattern of 
correlation . 
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CLAIMS 

1- A method of constructing a deterministic model 

of the diffusive type, which allows the multi- lithological 
filling of a sedimentary basin to be simulated, in which, 
5 on the basis of known field data pertaining to the 
architecture of the basin and measurement data as well as 
well logging data or seismic data, a set of input data is 
set up relating to an accommodation structure brought 
about by subsidence and eustasy, the inflow of fluviatile 
0 or marine sediment and the transport thereof, and physical 
parameters such as the diffusion coefficients of different 
lithologies, wherein, by a process of iteration, it 
consists in: dividing the basin into grids of regular 
dimensions; solving large-scale diffusive transport 
5 equations in accordance with an explicit pattern of finite 
volumes with constant time intervals, so as to simulate 
the flux of each lithology deposited on each grid; and 
comparing the results of the simulation with the field 
data and modifying the input data step by step by a 
process of inversion. 

2. A method as claimed in claim 1, wherein the 
inversion stage is performed by a process of trial and 
error . 

3. A method as claimed in claim 1 or 2 , wherein the 
large-scale diffusive transportation of each lithology is 
modelled by the relationship: 



for each lithology i 




diffusion => 


Q = v K E 

*~- 1 K i X V i water 




'5, 


flux of the lithology 


I (in m 2 /s) 




lithology I content (dimensionless) 


where * K, 


diffusion coefficient 


relative to the lithology I (in m'/sl 




water efficiency (dimensicnless) 


th 


altitude of the ground surface ( in m> . 




and K. = K (;yv-: 
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4 . A method as claimed in one of the preceding 
claims, wherein the littoral transport of each lithology 
is modelled on a large scale by the relationship: 



for each lithoicgy i 



littoral drift => 



Q y = K,..EC g f(b)*inacosa r 



Sb/bp 

where f (b) = (e * c S P 



if b<bji 
if b>bp 



4 E 



flux of the lithology i induced by littoral drift {in m : /s) 
efficiency of the transport of the lithology i by the waves 



density of the wave energy 



lc ? rate of propagation of the wave energy 

! f (b) flux distribution function as a function of the bathymetry 
(a angle of approach of the waves 



5. A method as claimed in one of the preceding 
claims, wherein the mult iple - 1 ithology sedimentary layers 
are defined as consisting of a superposition of 
lithologies with porosities regarded as being independent 
of one another. 

6 . A method as claimed in one of the previous 
claims, wherein the introduction and production of 
sediments is modelled using a relationship such as: 



for each lithology i 



(p. 



production rate of the lithology I (in m.<si 



maximum production rate of the litboiogy 1 tin m/sj 



I 



I 



I influence of B. 
! the bathymetry 



fo 

lb 

lb. 
le P 



if b<0 



if 0<b<b t 
if b>b 



IB. influence of the bath>-tnetry < dimens icn les s « 

lb oathymecry '.ir. irt 

lb 4 bathymetn:: threshold ':r. r.; 

1(3. attenuation coefficient, nr. n * ■ 



{ 

i influence of F. 

I the _1 it he logy 
I 



! 1 

le V; 



if Q.<S. ; 

s " : if C.>S. 



I F, , sensitivity of lithology i :c Iitr.oiogy 
flux : 

IQ. flux of tr.e litr.olcgy ; < vr. ra'.'s* 

; S, : inhibition stare threshold !:n t/ s. 
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7. A method substantially as hereinbefore described 
with reference to the figures. 
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